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PREFACE 

The work described in this report was conducted by the Radar and 
Optics Division of the Environmental Research Institute of Michigan 
(ERIM). The work was supported by the Jet Propulsion Laboratory 
(JPL) Contract No. 956105, a subcontract under JPL's prime NASA con- 
tract NAS7-1C0. The report summarizes the entire contract perfor- 
mance period of August 1981 to April 1982. The JPL technical monitor 
for this report was Dr. Daniel Held. 

This report discusses only Task II of the contract. Task I, 
which was concerned with a technical review of the JPL User's manual 
for the Seasat SAR Data, has been completed and sent to Dr. D. Held 
under separate cover. 

Co-principal investigators for this program were Drs. Robert 
Shuchman and David Lyzenga. Mr. James Marks contributed the section 
on Image Scale Factors for ERIM-generated SEASAT SAR imagery. Mr. 
Alex Klooster performed the analysis on image film radiometric mea- 
surements, and Dr. D. Lyzenga carried out the study on radiometric 
effects of Doppler spectrum variations. 
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EVALUATION OF OPTICALLY PROCESSED 
SEASAT SAR DATA 

1 

INTRODUCTION 

This informal information report describes the results of three 
experiments which were performed to evaluate some of the characteris- 
tics of ER IM optically processed Seasat SAR data. Two previously 
reported engineering studies (Shuchman, et al., 1978, 1981) have 
evaluated various aspects of the Seasat SAR system and its products, 
including more detailed analyses of the data collection and recording 
systems, the image quality, and system radiometric effects. 

The first experiment performed during the present study was in- 
tended to determine how accurately the scale of an optically proc- 
essed image could be predicted using only the data recorded on the 
signal film and the information provided on the auxiliary data list- 
ing for the Seasat orbit considered. A set of images was made for a 
portion of Rev. 651 with a calculated scale of 1:250,000 and the re- 
sults were compared with standard topographic maps to determine the 
actual scale of the images. 

The second study dealt with the practical considerations involved 
in making radiometric measurement: from Seasat image films. The 
feasibility of using calibrated point targets, i.e., corner reflec- 
tors, to measure the radar cross section of distributed scenes is 
evaluated and recommendations are made to improve the accuracy of 
th's procedure. 

The third study explored the radiometric effects of Doppler spec- 
trim variations due to earth rotation and spacecraft altitude 
changes. These effects are calculated and compared with measurements 
using data from Seasat Rev. 757, and possible correction procedures 
are discussed. The results of this study are especially relevant to 
the interpretation of data contained in long swaths of Seasat 
imagery. 
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IMAGE SCALE DETERMINATION 
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This study was a test of the accuracy with which the scale of a 
given Seasat image can be predicted using only the system parameters 
contained in the Sensor Data Record (SDR), the properties of the film 
recording system as inferred from an examination of the signal film, 
and the magnification settings of the optical processor. The study 
was carried out by setting up the optical processor so that output 
imagery with a nominal scale of 1:250,000 was produced and then com- 
paring this imagery with a 1:250,000 scale map. This section of the 
report describes the method used to determine the processor settings 
and the results of the scale comparisons. 

2.1 PROCESSING METHOD 

The data used in this study were collected during Seasat Rev. 
651 near Louisville, Kentucky. The data were read from High Density 
Digital Tape (HDDT) and recorded on signal film using the ERIM/Seasat 
film recording system at the Applied Physics Laboratory (Shuchman, 
et al., 1978). This recording system has a nominal sweep speed of 
500 m/sec and a film drive speed of 32.9 mm/sec; however, in order 
to determine these parameters more exactly, a set of measurements 
was made on the signal film itself. 

To determine the recorder sweep speed, the pattern resulting from 
a 10 MHz test signal recorded on the signal film was examined. By 
counting cycles with a toolmaker's microscope, the spacing was found 

5 

to be 5.40 x 10 m/cycle. Multiplying this by the frequency 
(10^ cycles/sec), a sweep speed of 540 m/sec was obtained. 

To determine the film drive speed, a calibrated grating with a 
line spacing nearly equal to that of the signal film was placed over 
the film, and the resulting Moire fringe pattern was observed. Using 
a grating with 49.98 cycles/mm, 8.5 cycles were observed over a 
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13.9 rrni length of the Moire fringe pattern. Thus, the spatial fre- 
quency of the signal film Is 

F $ - 49.98 - ~] 3 ■ 49.37 cycles/nm . (1) 


01 v 1 ding this Into the pulse repetition frequency of 1646.75 Hz, the 
film drive speed was determined to be 33.36 mm/sec. 


Once these recording system parameters are known, the scale fac- 
tors of the signal film can be determined with information from the 
sensor data record. The azimuth scale factor can be calculated from 
the film speed if the swath velocity is known. This velocity was 
obtained from the sensor data record for each quarter-swath, as shown 
ir. Table 1. Using these values, the along-track scale factor of the 
signal film is 201,926 for swath 1, 201,920 for swath 2, and 201,915 
for swath 3. Swath 4 was not processed during this experiment. 


The slant range scale factor of the signal film is obtained from 
the sweeo speed V $ using the equation 


c 


2.998 x 


*8 


10 w m/sec 
m/sec) 


277,600 


( 2 ) 


The ground range scale factor is equal to the above divided by the 
sine of the incidence angle. Thus, at the center of swath 1 the 
ground range scale factor of the signal film is 277,600/sin 20.698* » 
785,420. Similarly, the ground range scale factors for swaths 2 and 
3 are 726,000 and 676,480 respectively. Note that the ground range 
scale factor changes continuously acros; the swath so that at the 
edges of swath 1, for example, the scale factors are 819,630 and 
754,130 ( i . e . , about *4 percent different from the scale factor at 
th» center of the swath). 

The scale factor of the image film is determined in the range 
direction by the processor range magnification, and in the azimuth 
direction by the ratio of the input film speed to the output 
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TABLE 1 

SDR DATA FOR PORTION OF REV. 651 CONSIDERED 


Swath Swath Veloc Ity Incidence Angle (at Near Range) 


1 6736.253 m/sec 19.797* 

2 6736.061 m/sec 21.559* 

3 6735.870 m/sec 23.365* 

4 6735.680 m/sec 25.093* 
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film speed. The range magnification of the ERIM optical processor 
is 1.0044. Thus, to achieve unity aspect ratio at the center of 
swath 1, the film speed tracking ratio must be 


R 


t 


785,420/1.0044 

— ftT7&5 


3.873 


(3) 


For swath 2, the tracking ratio is 3.580 and for swath 3 the tracking 
ratio is 3.335. Image films ware generated for each of these swaths 
using these processor settings. 


Finally, the image films were printed using the appropriate mag- 
nifications in the enlarger to obtain a nominal scale of 1:250,000. 
The appropriate magnification is the ratio of the image film scale 
to the desired print scale, i.e., 3.128 for swath 1, 2.891 for swath 
2, ana 2.694 for swath 3. An example of the imagery for swath 1 is 
shown in Figure 1. 


2.2 EVALUATION 

The actual scale of the output prints was evaluated by comparing 
them to 1:250,000 ^cale U.S. Army Corps of Engineers map number 
NJ16-6. Features along the Ohio River, which runs in approximately 
the range direction were located near the edges and center of each 
swath. The acrcss-track scale factor was computed separately for 
each half of the swath by comparing the distances between features 
near the center and at the edges of the swath, and for the entire 
swath by comparing distances between features near opposite edges of 
the swath. The results are summarized in Table 2. Note that there 
is a consistent pattern of larger image scale factors for the left 
side of each swath than for the right side. This is in agreement 
with expectations since the ground range scale factor depends on the 
incidence angle and t^e processor was set up to yield unity aspect 
ratio at the center of each swath. 
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FIGURE 1. PORTION OF REV. 651 IMAGE USED FOR SCALE EVALUATION. 
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TABLE 2 

OBSERVED DIFFERENCES IN ACROSS-TRACK SCALE FACTORS 
BETWEEN SAR IMAGE AND MAP 


Swath 

Left Side (%) 

Right Side (%) 

Entire Swath (X) 

1 

+2.8 

-3.3 

-0.4 

2 

+0.8 

-6.1 

-2.7 

3 

-1.5 

-5.7 

-3.5 
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There also appears to be an overall decrease in the image scale 
factor from swath 1 to swath 3. The reason for this is not known 
but may be due to errors in the SDR incidence angle data* errors in 
processor settings, or inaccuracies In the scale factor measurements. 
The along-track scale factor was found to be accurate to within 0.3 
percent of the map scale factor. 

2.3 SUWARY 

This study was designed to test the accuracy with which a speci- 
fied image scale factor could be obtained in optically processed 
Seasat imagery. The results indicate that by using only the system 
parameters given in the SDR, the properties of the film recording 
system, and the magnification settings of the optical processor, an 
image scale factor can be predicted to within -*0.3 percent for the 
azimuth (along-track) direction. 

In the range ( across-track) direction, the image scale factor is 
observed to deviate by as much as 6 percent from the predicted value 
(see Table 2). A -±4 percent deviation was expected for the data 
studied, assuming that the specified scale factor was produced at 
the center of the swath. Such deviations can be expected in any 
optically processed imagery because the processor cannot compensate 
for the decrease in the ground range scale factor due to Incidence 
angle variations within each swath. The difference between the ob- 
served 6 percent deviation and the expected 4 percent deviation is 
apparently due to an overall decrease in the image seal"* factor from 
swath 1 to swath 3. This change may be due to errors in processor 
settinqs, interpretation of SDR incidence ancle data, or tc measure- 
ment inaccuracies. 
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The microwave radiometric characteristics of objects resolved by 
a Synthetic Aperture Radar (SAR) may be ooserved in the output image 
of the SAR. Radar cross section measurements are usually made photo- 
metrically in the image plane of the optical processor. However, it 
is sometimes necessary to make radiometric measurements on the output 
imagery as recorded on film, for example where the investigator does 
.not have access to an optical processor. This section of the report 
discusses the problems involved in extracting backscatter cross sec- 
tion values from SAR image films, and makes recommendations for im- 
proving this procedure. 

3.1 BACKGROUND 

The usual method of making cross section measurements in a SAR 
optical processor is to record the electrical current generated by a 
photomultiplier tube (PMT) intercepting the light flux passing 
through a suitable aperture in the image plane. When measuring the 
cross section of a point target, the aperture is chosen to pass only 
the light in the point image. Measurements of average reflectivity 
for distributed targets .are made by adjusting the aperture to pass 
the flux representing the area to be averaged. If the scale of the 
image plane is known, the aperture size in radar space may be deter- 
mined and measurements of reflectivity per unit area can be computed. 
Absolute measurements may be inferred if objects of known cross sec- 
tion, such as calibrated corner reflectors, are included in the 
image. 

When the image produced by the SAR optical processor is recorded 
on photographic film, the photometric properties of the image are 
transformed to density variations on the recording film. If the re- 
cording characteri sties of the image film are precisely known, it is 
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possible to calculate the original photometric properties of the 
image and thus the radiometric properties of the imaged area. These 
characteristics can be established in part by a step wedge exposure 
on the image recording film. The additional information needed to 
obtain absolute radar cross section measurement from the image film 
is the film density for an object or area of known radar cross sec- 
tion. The circulation of radar cross sections for distributed tar- 
gets is discussed in the following section. Additional considera- 
tions involved in the use of point targets, such as corner reflec- 
tors, for calibration are discussed in Section 3.3. 

3.2 RADAR CROSS SECTION MEASUREMENTS 

The relative radar reflectivity of two radar image areas may be 
defined as the ratio of the cross section per unit area (o Q ) for 
each location, i.e., 

Backscatter Ratio ■ o o2 la Q ^ (4) 

where o q2 is the reflectivity of the test area, and o Ql i$ 
reflectivity of a reference area. 

The reflectivity of the terrain is related to the image inten- 
sity, radar calibration constant, and collection geometry by 

o Q - IC cos 6 (5) 

where I is the image intensity, C is the radar calibration factor, 
and 6 is the depression angle. 

The calibration factor C includes the effect of antenna gain and 
range to the target as well as other factors. It therefore generally 
varies over the radar image and must be taken into account when ob- 
jects at different ranges and depression angles are compared. The 
cosine factor appears in the o Q expression because the SAR image 
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has coordinates corresponding to the slant range plane and o Q is 
defined in ground coordinates. 

The reflectivity ratio may then be related to the image intensi- 
ties as 



°o? j| ! 2 C 2 cos *2 

° 0 1 I-jC-, cos 0.| 


( 6 ) 


When the image is recorded on film, the intensities I-j and 
12 are converted to corresponding image densities D-| and D 2 by 
action of the film exposure E and the photographic development of 
the film. By definition 

E - It (7) 

where t is the effective exposure time of the recording process. 

The usual way of relating the exposure and density characteris- 
tics of a film is by the D vs. log E curve. An example is shown iri 
Figure 2. 

When the density and log exposure scales are equal, the slope of 
the straight line portion of the curve describes the contrast or 
"gamma" of the film; i.e., y - Tan a where a is the angle of the 
straight line portion with respect to the horizontal axis. When 
D-j and D 2 fall on the straight line portion of the D vs. log E 
curve, the image intensity ratio can be found by using the relation- 
ship developed below: 

0, - 0 2 =. T [log E, - log E 2 ] (8) 

. i log (E,/E 2 ) 

- t log (Ijt/Ijt) 

= y log (I,/I 2 ). 


Then 
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FIGURE 2. TYPICAL DENSITY VS. EXPOSURE PLOT FOR PHOTOGRAPHIC FILM. 



FIGURE 3. EFFECTIVE FILM "GAMMA' 1 FOR POINTS NOT ON THE LINEAR PART OF 
THE CURVE. 
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( D 1 - D 2 )/ y - log (X 1 /I 2 ). (9) 


This intensity ratio may then be related to the relative back- 
scatter ratio o^/o^ given in Eq. (6). The ratio in dB is given 
by 


^2 cos ^2 


10 log -2i - 10 log r r ■ - c - i - . 
°ol ’l C l C0S *1 


( 10 ) 


Substituting for the intensity ratio, the density relationship in 
Eq. (9), the backscatter ratio in dB is given by 

o 9 D ? - D, C ? cos 

10 log — - 10 log + 10 log -J- + 10 log — A . (11) 

a i y cos v i 


For situations where the depression angle t and the calibration con- 
stant C are essentially the same for the two areas to be compared, 
the backscatter ratio in dB is approximately 


10 log — = 10 log 
°ol 



( 12 ) 


If the densities D 2 and are not on the straight line por- 
tion of the curve, the intensity ratio may still be found using the 
above relations if y is computed as the slope of the line joining 
points D 2 and D ( on the 0 vs. log E curve. This is shown graph- 
ically in Figure 3. 

When the resulting value of y is significantly different from 
that of the straight line portion of the 0 vs. log E curve, the accu- 
racy of the ratio measurement will be reduced. For this reason the 
film recording exposure should be chosen to place the densities of 
all the areas of interest on the straight line portion of the D vs. 
log E curve if possible. 

The absolute value for the cross section of a test area, 0 ^, 
can he obtained only by making measurements relative to a reference 


\ 
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area with known radar cross section. Frequently, point targets such 
as corner reflectors are used for calibration. However, the measure- 
ment of peak density for such point objects on the image film Is 
difficult because of the small size of the images and the lack of 
suitable instrumentation for making point density measurements. 
Microdensitometer scans may be used, but multiple scans must be made 
to be sure that the peak density is recorded. The densities produced 
by point image exposures may also have sensltometrlc characteristics 
that differ from those made from large area exposures such as pro- 
duced by a standard step table. A measurement of this effect is 
described in Section 3.3 In addition, the reflector responses will 
often be saturated in images where the exposure level is chosen to 
record the average terrain at normal density levels. Because of the 
factors cited above it is extremely difficult if not impossible to 
calculate accurately the relative energy content of reflector re- 
sponses by measuring their corresponding image densities. Fortun- 
ately, the measurement of average exposure levels for distributed 
radar returns usually relates directly to standard photographic 
sensltometer exposures, as shown in Section 3.4, and measurement of 
relative average backscatter may be readily made using standard 
photographic instrumentation. Absolute measurements of average back- 
scatter are then possible if a reflectivity standard for the radar 
image, established by the usual photometric means in the optical 
processor, is recorded as part of the imagery package. 

3.3 IMAGE POINT DENSITY MEASUREMENTS 

Point density vs. exposure level measurements were made for 
Seasat corner reflector images of the Goldstone reflector array as 
recorded on signal film No. 286 of Rev. 882. The measurements were 
made using imagery recorded cn Kodak 3414 film and a modified form 
of the Joyce-Lobel microdensitometer. The radar resolution was 
approximately 10 meters in azimuth by 25 meters in ground range. 
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The image scale was approximately 700,000:1 yielding an image film 
resolution of about 14 ym by 35 urn. It was not possible to make 
accurate peak density measurements of these small spots using the 
standard optical arrangement of the Joyce-Lobel microdensitometer 
because it was impossible to observe the alignment of the spot image 
at the aperture screen. To circumvent the problem, the image was 
projected with approximately a 10— X magnification into a 50 ym Gamma 
Scientific photometer probe by means of a mirror, using the illumina- 
tion source and scanning table of the Joyce-Lobel instrument. The 
scanning table of the microdensitometer was used to position the film 
so that the point images were properly aligned with the 50 ym probe. 
Densities (D) were then calculated by recording the light level indi- 
cated on the Gamma Scientific autophotometer and then using the 
relation 

D - log (B/P) (13) 

where B is the background reading with film removed, and P is the 
minimum light reading through the reflector image. 

The point density vs. exposure characteristic was established by 
generating a sequence of images on the ERIM precision optical proc- 
essor using 3 dB light level increments. A standard sensitometer 
exposure was included to establish the D vs. E characteristic for 
the film. 

The light level measurements used to compute the image densities 
are shown in Table 3. The measurements were made difficult by the 
presence of a granular "anti-stick" coating on the back of the film. 
The random orientation of the relatively large grains lead to a dif- 
ferential refraction of the light beam away from the measurement 
path. Consistent density measurements were possible only after the 
back side of the 3414 film was covered with immersion oil and a glass 
plate. 
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TABLE 3 

CORNER REFLECTOR IMAGE DENSITY MEASUREMENTS 
Background Measurement B = 20 (In Same Units As P, Below) 



Reflector 

1 

Reflector 

2 

Exposure 

P 

D 

P 

D 

0 

52.0 

0.32 

98.0 

0.09 

+3 dB 

20.0 

0.78 

45.0 

0.43 

+6 dB 

4.4 

1.44 

14.0 

0.93 

+9 dB 

0.96 

2.10 

2.3 

1.71 

+12 dB 

0.42 

2.46 

0.58 

2.32 

+15 dB 

0.32 

2.57 

0.36 

2.52 
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The point densities calculated above are plotted in D vs. log E 
format In Figure 4. The standard sensitometer densities as measured 
with the McBeth densitometer are plotted on the same graph. The 
slope of the D vs. log E curve for the two point image responses is 
approximately 2.2 as shown on the graph. This is significantly 
hiaher than the 1.55 slope of the standard sensitometry curve. Addi- 
tional step wedge density measurements were made using the point 
density instrumentation. The slope of the resulting curve as shown 
in Figure 4 was slightly higher than that for the McBeth measure- 
ments, but still much lower than that of the point measurements. 

The reason for this difference is not known, but may be due to 
adjacency effects in the film development process. These effects 
are caused by differing rates of depletion of the developer for point 
images and diffuse images. 

3.3 DISTRIBUTED REFLECTANCE MEASUREMENTS 

In this part of the study, measurements were made cf the image 
film density for a uniformly reflecting area as a function of the 
exposure level. This was done to confirm that the sensi tometri c 
properties of the film were the same for such distributed reflectors 
as for the standard sensitometer exposures. A simulated signal film, 
consisting of a Mylar diffuser, was used in the optical processor. 
The proc^sso** bandwidths were set up for Seasat process'^ and an 
image film was exposed using various laser illumination levels. 
Prior to development of this *ilm, a standard set of sensitometer 
exposures was also added to the film. 

After development, the density of the sensitometer step wedges 
was measured and plotted versus the exposure. The density of the 
simulated images was also measured and plotted versus the laser illu- 
mination level. The resulting plots, shown in Figure 5, indicated 
similar characteristics (i.e., y) for the step wedges and the test 
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exposures. Thus it appears to be valid to apply standard sensito- 
metric techniques to distributed targets in SAR imagery. 

3.4 SUMMARY AND RECOMMENDATIONS 

The conc 1 usicn which follows from these measurements is that it 
is not feasible to use point targets, such as corner reflectors, to 
calibrate the image film for distributed targets since the film re- 
sponse is different for point targets and distributed targets. One 
possible solution of this problem is to make measurements of the 
point target response in the optical processor and use these measure- 
ments to infer the c-oss section per unit area, o Q , of a jniform 
area in the scene. Since distributed targets do appear to have the 
same densitometric properties as the standard sensitometer exposures, 
a uniform area of known a Q could then serve as a secondary calibra- 
tion reference for the image film. 

The accuracy of both relative and absolute backscatter measure- 
ments is generally reduced as the time between data collections for 
the two areas to be compared increases. Even during a single data 
pass the radar calibration factor C (Eq. (5)) may vary for identical 
radar ranges as a result of combinations of the following: 

1. Receiver gain changes 

2. Transmitter power variations 

3. Signal film density variations 

4. Azimuth (Doppler) spectrum wander due to antenna orientation 
changes and electronic offset errors 

5. Antenna gain variations due to radar platform roll. 

If the effects of these variations are net properly accounted for, 
serious backscatter measurement errors will result. In particular, 
failure of the raoar motion compensation system to adequately control 
azimuth spectrum wander often produces "banded" imagery with gross 
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changes of image intensity over relatively short along-track dis- 
tances. Seasat radar data are also affected by Doppler spectrum 
shifts due to earth rotation effects as described in the following 
chapter of this report. Accurate backscatter calculations may be 
made under the above conditions if their effect on the value of the 
calibration constant are precisely known. If they are unknown, back- 
scatter measurements will be equally uncertain. 


Terih 


RADAR AND OPTICS DIVISION 


4 

RADIOMETRIC EFFECTS OF DOPPLER SPECTRUM VARIATIONS 

The along-track radiometric stability of the Seasat SAR is an 
issue of concern to those attempting to correlate the SAR image in- 
tensity with geophysical parameters, such as winds, over large areas. 
The study reported in this section deals with one aspect of this 
problem, namely the effect of Doppler spectrum variations on the 
radiometric properties of the SAR image. The discussion is in the 
context of optically processed imagery since the ability of optical 
methods to process long swaths in a continuous manner makes these 
methods attractive for this application. The effects of Doppler 
spectrum "wandering" on digitally processed data must also be consid- 
ered when comparing successive frames of imagery, but these effects 
are beyond the scope of the present study. 

4.1 BACKGROUND 

The azimuth, or Doppler, spectrum of the signals received by a 
SAR is dependent on the antenna gain pattern, the antenna look direc- 
tion, the platform velocity, and the radial velocity of the objects 
in the scene. For Seasat, the received Doppler spectrum has a width 
of about 900 Hz. Because this spectrum is, in effect, sampled at 
the pulse repetition frequency (PRF) the actual spectrum is not dis- 
tinguishable from a set of "alias" spectra separated in frequency by 
multiples of the PRF. Approximating the antenna gain pattern by a 
Gaussian function, the received signal spectrum can be written as 



n=-oo 


where o is the radar cross section of the surface, B = 900 Hz is 
o 

the 3 dB bandwidth of the Doppler spectrum, f c is the peak fre- 
quency, and f = 1647 Hz is the pulse repetition frequency ( p PF). 


25 


PRECEDING PAGE BLANK NOT FILMED 


RADAR AND OPTICS DIVISION 


When this signal is recorded on film, the higher frequencies are 
attenuated, or recorded with less efficiency than the lower fre- 
quencies. This effect is described by the modulation transfer func- 
tion (MTF) of the recording system. Measurements in the optical 
processor indicate that the shape of the MTF can be approximated by 
the sum of a Gaussian and an exponential function, i.e.. 



where B f> ■ 600 Hz Is the receding system bandwidth. Thus, the spectrum 
of the signals recorded on the signal film (as measured in the optical 
processor) is given by 

S Q (f) - M(f)S.(f) ♦ N(f) (16) 

where N(f) is the noise spectrum, which has approximately the same 
shape as the MTF. Examples of measured spectra for Seasat Rev. 757 
are shown in Figures 6-10. Note that the effects of the MTF are 
clearly visible in these measurements, causing the alias spectra to 
be of unequal amplitude and the peak frequencies to be separated by 
less than the PRF. Note also that the peak frequencies vary as a 
function of the distance along the swath during this pass. The im- 
plications of this spectrum shift for the radiometric calibration of 
the data will be discussed in Section 4.3. 

4.2 DOPPLER SPECTRUM MEASUREMENTS 

In order to measure the spectral shifts during Rev. 757, the 
spectra shown in Figures 6-10 were divided by the MTF. Using the 
measured MTF, the data appeared to be over-corrected, as indicated 
by the amplitudes and separation of the real and alias spectra. 
Using a modified MTF of the same form as Eq. (15) but with B r , 
900 Hz, a better correction was obtained. The corrected spectra 
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FIGURE 6. AREA A 


FIGURE 7. AREA B 


FIGURE 8. AREA C 
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using this MTF are shown In Figures 11-15, and the peak frequencies 
of the real and alias spectra are shown in Table 

Shifts in the peak frequency of the Doppler spectrum along the 
swath are due to earth rotation effects and spacecraft attitude vari- 
ations. The Doppler shift due to the earth's rotation is given by 

2 SIR. 

tf — cos (L) sin i sin e (17) 

v A 

where ft is the earth's angular velocity (7.272 x 10“^ rad/s), R g 
is the earth's radius (“6370 kmd), x is the radar wavelength (23.5 
cm), L is the latitude, t is tne radar look direction (from North), 
and e is the incidence angle (-20*). 

The Doppler shift due to spacecraft attitude variations is given 
approximately by 

af c ■ -f (cos e sin P - sin e sin Y) (18) 

where V is the swath velocity (nominally 6800 m/s), P is the Ditch 
angle, and Y is the yaw angle. P is positive for nose-up rotations 
of the spacecraft and Y is positive for clockwise rotations (looking 
down). These angles are given in the auxiliary data listing for the 
Seasat pass under consideration. 

The latitude variation and changes in pitch and yaw for a portion 
of Rev. 757 are shown in Figure 16. Also shown in this figure are 
the calculated Doppler shifts due to the earth's rotation and the 
changes in spacecraft attitude along the pass. The total Doppler 
shift is shown in Figure 17, along with the peak frequencies obtained 
from the optical processor measurements. Note that the trend is 
predicted quite well by the above equations. There is an offset of 
about 200 Hz between the calculated and observed Doppler frequencies, 
but this is to be expected since there are Doppler shifts introduced 
during the transmission and recording of the data which are not 
accounted for in these calculations. 
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input signpl - 20 m from end i»*>ut sicnpl -15 m from end 




FIGURE 11. CORRECTED DOPPLER FIGURE 12. CORRECTED DOPPLER 

SPECTRUM FOR AREA SPECTRUM FOR AREA 

A, REV. 757. B, REV. 757. 


input siGNBi. -10 m from end 



FREQUENCY (Hz) 


FIGURE 13. CORRECTED DOPPLER 
SPECTRUM FOR AREA 
C, REV. 757. 
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input sicw -5m from end 



FIGURE 14. CORRECTED DOPPLER 
SPECTRUM FOR AREA 
D, REV. 757. 


input sicw.- end of pass 



FIGURE 15. CORRECTED DOPPLER 
SPECTRUM FOR AREA 
E, REV. 757. 
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TABLE 4 

DOPPLER PEAK FREQUENCIES MEASURED FOR 

REV. 757 

Location 

T ime 

f (actual spectrum) 
c 

f (alias spectrum) 

A 

22:37:36 

400 Hz 

-1300 Hz 

B 

22:40:13 

700 Hz 

-950 Hz 

C 

22:42:49 

900 Hz 

-730 Hz 

D 

22:45:25 

1050 Hz 

-500 Hz 

E 

22:48:01 

1250 Hz 

-400 Hz 


\ 
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1400 J 


1200 J 



TIME (MIN) 


FIGURE 17. CALCULATED DOPPLER SHIFTS (SOLID LINE) 
DUE TO EARTH ROTATION AND SPACECRAFT 
ATTITUDE CHANGES, AND OBSERVED PEAK 
FREQUENCIES (CIRCLES) FOR SEASAT REV. 
757. 
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4.3 RADIOMETRIC EFFECTS 

The effects of changes in the Doppler frequency on the radio- 
metric calibration cf the data depend on how the data are processed. 
Normally, an aperture is placed in the frequency plane of the proc- 
essor to block out the alias spectrum, in order to avoid azimuth am- 
biguities in the image. If this is done, the effective gain of the 
system will vary significantly along the swath, even if the location 
of the aperture is moved so as to track the changes in the spectrum 
location. This effective gain, defined as the ratio of the output 
power to the input power within the frequency aperture, is plotted 
versus the peak Doppler frequency in Figure 18 for an aperture width 
of 900 Hz centered about the peak frequency. The same gain factor 
is plotted versus time for Rev. 757 in Figure 17. It is apparent 
that this is a very significant effect which should be accounted for 
if the image intensity, or image film density, is used to infer the 
radar cross section of the surface* over a long swath. 

The radiometric effects of Doppler shifts are much less pro- 
nounced if measurements are made in the frequency plane of the proc- 
essor using a wide enough aperture to pass both the primary and one 
alias spectrum, since a decrease in the amplitude of one peak is 
accompanied by an increase in the amplitude of the other. Note that 
^ hi s method effectively bypasses the image formation process and thus 
results in a much lower spatial resolution. There is still a slight 
variation in the effective gain of the system using this method, as 
shown in Figures 20 and 21, but the variation is small compared to 
other sources of error. 

An alternative to this procedure would be to make explicit cor- 
rections for Doppler spectrum variations. This could be done by 
measuring the peak Doppler frequency at several points along the pass 
and calculating the gain as was done here. Such a correction would 
be especially important for measurements made on conventionally pro- 
cessed image films over large swath lengths. 
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FIGURE 18. RATIO OF OUTPUT TO INPUT POWER WITHIN A 
FREQUENCY INTERVAL OF 900 HZ CENTERED AT 
THE DOPPLER PEAK FREQUENCY. 
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FIGURE 19. RATIO OF OUTPUT TO INPUT POWER FOR A 

SEGMENT OF REV. 757, USING A FREQUENCY- 
PLANE APERTURE OF 900HZ WIDTH CENTERED 
AT THE ACTUAL DOPPLER PEAK FREQUENCY. 
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FIGURE 20. RATIO OF OUTPUT TO INPUT POWER FOR A 
FREQUENCY INTERVAL FROM -PRF TO +PRF 
VERSUS DOPPLER PEAK FREQUENCY. 
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FIGURE 21. RATIO OF OUTPUT TO INPUT POWER FOR A 

SEGMENT OF REV. 757, USING A FREQUENCY- 
PLANE APERTURE FROM -PRF TO +PRF. 
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4.4 SUMMARY 

Using optical processing techniques, variations in the Doppler 
spectrin of the received signal are observed to significantly affect 
the radiometric characteristics of Seasat SAR image data. These 
effects are sufficiently important that they should be accounted for 
if the image intensity is to be used to infer radar cross section 
over a long image swath, as, for example, in the measurement of sur- 
face wind speeds. It is found that these effects are much less pro- 
nounced if measurements are made in the frequency plane of the proc- 
essor using an aperture which is wide enough to pass both the primary 
and one alias sp'-'.*rum. These effects can also be corrected by 
determining a rali. region curve of the effective gain for the pass 
from measurements of the peak DoDpler frequency as a function of 
time. This correction method is especially applicable to intensity 
measurements made on conventionally processed Image films over large 
swath lengths. Determination of proper MTF is important in the lat- 
ter correction method, but is not needed in the frequency plane 
method . 
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